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ABSTRACT 


In this thesis, large deflection plate compatibility 
and equilibrium differential equations are derived using finite 
deformation theory including the non-linear five constant stress- 
strain relationship. In addition, the assumptions of classical 
large deflection plate problems under classical elasticity 
theory are discussed in connection with the more general 


approach proposed by this thesis, 
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INTRODUCTION 


The approach to the large deflection plate problem pro- 
posed by this thesis, as suggested by the title, is one utilizing 
non-linear elasticity theory, or theory dealing with finite defor- 
mation of elastic solids. In reference |1] Francis D. Murnaghan 
provides an excellent discussion of the theory, and derives the 
basic stress-strain relationship used in this thesis (see appendix 
B). The use of this basic non-linear stress-strain relationship, 
sometimes called Murnaghan's law |2], places the scope of this 
thesis more appropriately into the field of rheology rather than 
structures per se (see Figure 1). 

In actuality, predicting how materials will behave in 
response to forces lies in the domain of rheology - literally 
the study of flow. The Society of Rheology uses the words of 
Heraclitus, M%oavT@ Cer (everything flows), as its motto, and 
books on rheology quote this philosophical conviction in their 
introductory pages. [3]. Because rheology deals with the flow 
of matter under the action of forces, by an obvious extension 


it also includes all deformations of materials by forces. Thus, 


l1] References are listed beginning on page 66. 


L2] Novozhilov states "The elastic law corresponding to it (the 
five constant stress-strain theory) is ordinarily called 
Murnaghan's law, although it was actually first proposed 
much earlier by Voigt in 1893. The first attempt to examine 
the stress-strain relation in a form different from Hooke's 
law was made by Bulffinger in a paper published in the works 
of the Russian Academy of Sciances in 1729" pp. 127. 
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as shown in Figure 1, rheological theories have as their model 

the classical theories of elasticity and hydrodynamics in which 

it is assumed that a material, whether solid or liquid, deforms 
linearly in response to a force. From the classical models have 
grown more and more general theories, first by abandoning classical 
distinctions between a solid and a liquid, and then by abandoning 
the requirement that the response be linear. The most general 
theory, as shown in Figure (1), is the theory of nonlinear visco- 
elasticity which includes all other theories as special cases. 

A theory of non linear visco-inelasticity exists (shown in Figure 


1 dashed) but is of historical interest only [3]. 
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FIGURE 1 


(From "Deformation and Flow" by Charles J. Lynch as appeared in 
International Science and Technology January 1966) 
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Since the end of World War II there has been a consid- 
erable change in the approach to rheology and to the formulation 
of rheological theories. The origin of this change is two-fold. 
First, certain flow phenomena, observed first in flame-thrower 
fuels and later in a wide variety of other fuels, were in direct 
opposition to classical hydrodynamic theory. The other factor, 
which was perhaps of greater significance in its contribution 
to the changing approach to rheology, at least to those primarily 
interested in structures, was the development of a continuum 
theory for elastic materials. 

Although the theory of finite deformation of an elastic 
solid applies equally to all elastic materials, most of the 
research accomplished to date has been in the rubber industry. 
The theory as applied to metallic structures has only recently 
become of practical interest with the increased emphasis on 
limit design, and the more stringent requirements of high per- 
formance aircraft and missiles, space craft, and deep submergence 
pressure hulls. In the majority of experimental research con- 
ducted so far, rubber is generally used because it has a large 
elastic range and its nonlinearities are therefore more easily 
measured. 

Vulcanized rubber differs from other elastic materials in 
the extent of its elastic deformation. <A rubber rod or strip 
may be stretched to four or five times its initial length without 
suffering permanent deformation. When rubber is subjected to 


Shear deformation and the corresponding shearing force is measured, 
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the relation between the shear strain and the shearing force 
is found to be approximately linear, just as with most elastic 
solids. On the other hand, when rubber is stretched by tensile 
forces, the relation between the tensile force and the elonga- 
tion shows considerable departure from linear behavior. 

While this may at first sight appear surprising, it can 
in fact be demonstrated mathematically that a material having 
a linear shear relationship cannot possibly have a linear tensile 
relationship. From this we can conclude that a true Hookean 
solid is only a convenient fiction - an elastic material cannot 
possibly have a linear relationship in both tension and shear as 
assumed in classical elasticity theory. The classical theory is 
only valid because the departure from linearity does not show up 
when the deformations are small - even a curved line appears 
straight if one is only concerned with a little piece of it. [3]. 

In finite elasticity theory, the material properties are 
characterized by an expression for the energy of deformation per 
unit initial volume, (p(n) ), L1] (see appendix B). This strain 
energy function and the associated material constants have been 
determined by experiment for various vulcanized rubbers and the 
actual manner in which stored energy devends on the deformation 
has been determined from them. Using this experimentally deter- 
mined stored energy, the results of other experiments with vulcan- 
ized rubber have been predicted with considerable accuracy, pro- 
viding an excellent verification of the theory ue Although 


the material constants of interest to the structural engineer 
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(modulus of elasticity and shear modulus plus the three addi- 
tional third order elastic constants shown in appendix B) have 
been evaluated for several structural materials as shown in 
Table I, the numerical values show a wide variation between 
investigators with a wide spread for some materials and much 
more work is required. 
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Finite elasticity theory has successfully predicted a 
number of effects that are not predicted by classical theory. 
For example, experimental observations of a rod of circular 
cross section subjected to simple torsion made by Polynting 
in 1909 with steel wires and in 1913 with vulcanized rubber 
rod showed that the rod will not only twist as predicted by 
classical theory, but will also elongate. As is often the 
case when knowledge available in the literature remains unused 
for some time, these results seem to have been forgotten until 
recently when the application of finite deformation theory 
predicted this result [3]. 

That Murnaghan's finite deformation theory |1] should 
be applied to plates and shells in the large deflection region 
has been proposed verbally by S. F. Borg for some time. In refer- 
ence [4], Borg, Hoppe and Kopchinski, and in references [5], [6], 
[7] and [8] Borg discusses large deflection of plates, including 
the derivation of plate compatibility equations from a finite 
deformation theory approach. This thesis is a continuation and 


extension of Borg's work as applied to plates. 
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GENERAL DISCUSSION OF NON LINEAR ELASTICITY 


Non linearity is introduced into the theory of elasticity 
dealing with isotropic materials in three ways. It can be intro- 
duced through the generalized strain tensor or strain displacement 
relation (appendix A), through the stress-strain relations (appen- 
dix B) or through the equations of equilibrium of a volume element 
of the body (appendix C). 

In the equations of equilibrium and the strain tensor 
relations, the retention of non-linear terms is conditioned by 
geometric considerations, or the necessity of considering the 
angles of rotation in determining dimensional changes of line 
elements and in formulating the conditions of equilibrium of a 
volume element. However, in the stress-strain relation, non- 
finear terms appear if the strain exceeds in magnitude, physical 
constants characteristic of the material and referred to as the 
limits of proportionality. It should be noted that the above 
discussion does not presuppose a linear elasticity origin. 
Actually, linear theory is developed by successively ignoring 
those terms or quantities which give non-linear characteristics 
to the resultant equations. For many materials, e.g. aluminun, 
the limit of proportionality can be quite low or non existant, 
and a finite deformation, continuum, theory is required where 
deformations are sufficiently large as to invalidate linear 
elasticity theory (Hooke's law). 

Thus, two types of non-linearity must be considered in 


large deflection problems, geometrical and physical. Since, in 
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general, smallness of angles of rotation of the body as a whole 
does not imply the smallness of elongation or shear strain of an 
infinitesimal element of volume of the body (and conversely), the 
geometrical and physical non-linearity can be regarded as inde- 
pendent of each other. This independence is useful for purposes 
of classifying problem types and for purposes of making reasonable 
assumptions to simplify the problem to enable an engineering 
solution. 

As a result, elasticity problems can be classed as one of 


four basic types [2]: 


(1) Elasticity Problems having both physical and geometrical 
linearity. 


In problems of this type, the effect of the angles of 
rotation are of the same order of magnitude as the elongations 
and shear strains, while the elongations do not exceed the limit 
of proportionality of the material. An example of this type of 
problem is the simple tension test when the stresses are main- 


tained below the proportional limit. 


(2) Elasticity Problems which are physically non linear 


but geometrically linear. 


For these problems, angles of rotation can be neglected 
in projecting the forces which act on a volume element and in 
determining strains. However,’ the elongations exceed the pro- 


portional limit and require a non-linear stress-strain relation. 
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The simple tension test becomes a problem of this type when the 


stresses in the rod exceed the proportional limit. 


(3) Elasticity Problems linear physically, but non- 


linear geometrically. 


In problems of this type, the angles of rotation are 
essentially large but the strains do not exceed the limit of 
proportionality. An example of this type is the bending of a 
thin steel strip, or the buckling of a slender column within 
the elastic range where the original shape and position is re- 


gained after removal of the load, i.e. no permanent deformation. 


(4) Elasticity Problems non linear both physically and 
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In problems of the fourth type, the strains exceed the 
limit of proportionality and the angles of rotation are so large 
that it is necessary to retain non-linear terms both in the stress- 
strain relation, and in the equations of equilibrium. Plate and 
Shell problems become examples of this type if the deformation 


is large and the stresses exceed the limit of proportionality. 





LISCUSSION OF PLATE PROBLEM 


A plate may be defined as a three-dimensional deformable 
body whose thickness is small compared to its least lateral 
dimension. Because the thickness is small compared to the 
other linear dimensions, classical elasticity theory treats a 
relatively thin plate in an essentially two dimensional manner. 

The linear small-deflection theory of plates, developed 


by Lagrange (1811) is based on the following assumptions |9]: 


1. points which lie on a normal to the mid-plane of the 


undeflected plate lie on a normal to the mid-plane of the deflected 
plate ($4 == 0 ys 


2. the stresses normal to the mid-plane of the plate, 
arising from the applied loading, are negligible in comparison 


with the stresses in the plane of the plate ( 62<< Sx , by ) : 


3. the slope of the deflected plate in any direction is 
small so that its square may be neglected in comparison with 
nw 


unity (curvatures YT, = Geane ; he ~ 3y* 


4, the mid-plane of the plate is a neutral plane: l.e. 
any mid-plane stresses or membrane stresses arising from the 
deflection of the plate into a non-developable surface may be 
ignored. 

In classical large deflection theory where deflection is 
of the order of magnitude of the plate thickness, and greater, 


the first three assumptions above for small deflection theory are 
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retained. However, the fourth assumption is no longer valid with 
large deflection, and middle surface stresses must be considered. 
As a result of the large deflection, the problem becomes 
non linear. The question is, does the problem become one of 
geometric non linearity as assumed in classical theory, or is it 
one of both physical and geometrical non-linearity as proposed 
by this thesis? Physical non linearity is currently introduced 
in the inelastic buckling problem as a variation in Poisson's 
ratio and variation of the Elastic Modulus in the so-called tan- 
gent/secant modulus theory. Gerard, [10] provides a good "state 
of the art" summary of this theory. However, it is highly possible 
that the correct approach lies in the determination of the addi- 
tional elastic constants introduced by finite deformation theory. 
Assuming, for the moment, that the large deflection plate 
problem is physically linear, but geometrically non linear, there 
1S some error introduced by the classical differential equations 
formulated by vonKkirman [11] due to his treatment of the problem 
in a two dimensional manner. This error is negligible for very 
thin plates but becomes significant for plates of finite thickness. 
To illustrate this error, the classical large deflection assumptions 
will be applied to the generalized large deflection strain tensor 
(appendix A). 
The assumption that points which lie on a normal to the 
mid-plane of the undeflected plate lie on a normal to the mid- 
plane of the deflected plate require that 2 ZS, P(e) oe 
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Thus the large deflection strain tensor derived in appendix A 


reduces to: 
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are neglected, which seems reasonable, but even without the 
above terms, there are terms remaining in the expressions for 
Vy, Yyz and @3 which are still significant. 
If we assume that the elements of the generalized strain 


tensor, as applied to a plate of finite thickness, are of the 


order of magnitude shown below, i.e. 


(1) Terms y and gw are of the order € and are 
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of the same order of magnitude as the derivatives of the 


shear forces Q, and Qy 
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If the Poisson's effect across the thickness of the plate is 

ignored, and oy is taken to be zero, the above strain tensor 

becomes identical to that obtained by Borg, Hoppe, and Kopchinski |4]. 
In the differential equations formulated by von Karman Lael 


the strain elements are assumed to be 
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By comparing the von Kdérman elements with the elements of the 

three dimensional plate strain tensor, it is shown that von Karman 
retains terms of the order a and yet neglects terms of the 

order € which are introduced by the shear strain elements Yxz2 
and Wyz , and also ignores the Poisson's effect across the 
thickness of the plate Qu » Thus, even with the assumption 
that the problem is linear physically and non-linear geometrically, 
a feeling for the magnitude of the error resulting from the two- 
dimensional approach can be achieved. This points out a significant 
value of the continuum mechanics approach to the problem. 

At this point, it should be noted thet the large deflection 
strain tensor was developed in terms of a Lagrangian curvilinear 
coordinate system which at all points remains parallel to the de- 
formed plate (see appendix A). The treatment of a developable 
plate problem in the Lagrangian coordinate system lends itself 
readily to inextensional bending theory, but becomes a major diffi- 
culty in dealing with a flat plate with extensional plate theory. 
Thus, the use of the plate strain tensor derived above assumes 
that the deflections are small enough that an Kulerian coordinate 
system can be used. Even though it is recognized that error 1s 
introduced by the assumption of an Eulerian coordinate system, the 
plate strain tensor does provide an improvement over the von Karman 
two dimensional approach as discussed above, and should lead to 
less error for any practical plate problem. Indeed, if point 
strains are sufficiently large to invalidate the above assumptions, 


it is probable that plastic deformation would also invalidate 
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the assumption of an elastic body and a non-isotropic or plastic 
theory would then become necessary. 

In reference [12] Bleich gives an excellent historical 
sketch of the evolution of the von Karman equations and subse- 
quent strain energy approaches to the large deflection plate 
problem by Timoshenko, Marguerre and Trefftz. Ina series of 
papers by Levy 113], Praele PS |e 116] and Eivae solutions of a 
theoretically exact nature were given to the von Karman differen- 
tial equation [11]. In Levy's solutions, deflections and normal 
pressure were expressed in the form of Fourier series and sol- 
utions were obtained for various cases of loading and support. 
Although the method 1s involved and laborious, it has been recog- 
nized aS one of the more accurate solution techniques available 
under the von Karman theory. However, the only plate solution 
known to the author using finite deformation theory incorporating 
the associated third-order elastic constants is a preliminary or 
exploratory study of a circular plate by Borg in reference wer 
Although the solution in reference [7] is close to Timoshenko's 
solution of a similar problem contained in reference L138], the 
approach to the two solutions was vastly different. Rather than 
using the classical approach of Timoshenko, Borg's solution is 
obtained from the non linear stress-strain relation developed 
in appendix B, This would tend to lend evidence to the validity 
of the Finite Deformation Theory approach to large deflection 
plate problems, particularly where both physical and geometrical 


non linearity are required, 
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Another non linear large deflection thin plate approach 
using the non linear inextensional bending theory is described 
by Borg in references [5] and [8]. In this approach, a theory 
analogous to the Murnaghan non linear elasticity theory is 
developed. In the non linear thin plate inextensional bending 
theory, it is assumed, essentially,,.that the large deflection 
is one which occurs without the development of membrane stresses 


and a large deflection form of the bending deformation relation 


is obtained. [8]. The non linearity occurs in the second derivative 


of the deflection terms and arises as a consequence of the analo- 
gous behavior of a thin plate to that of an elastic body in gen- 
eral, and in a sense, bears the same relation to linear thin 
plate theory that the elastica (column) solution bears to the 
Kuler column theory. 

If the deflection of a plate is not small, the assumption 
regarding the inextensibility of the middle surface of the plate 
holds only if the deflection surface is a developable surface, 
therefore, there is some question about the generality of the 
inextensional theory. Although the non linear moment-—curvature 
tensor equations analogous to the non linear stress-strain rela- 
tions, as set down by Borg [8], provides a possible solution to 
the non linear plate problem, it would seem that the most general 
differential equations would result from the non linear stress- 
strain relation, where the classical equations of equilibrium 
incorporate non linear stresses, and the Compatibility equation 


results from the non linear stress-strain tensor equation, In 
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other words, the large deflection plate problem should be 

treated as one which is non linear both physically and geometri- 
cally, and ‘incorporates the third order elastic constants. 

Borg [6] has shown that the order of magnitude of the three 
additional elastic constants resulting from Murnaghan's stress- 
strain relation are of the same order of magnitude as the elastic 
constants G & E (see Table 1, page 5). As a result, their 
inclusion in the large deflection plate problem may be signi- 
ficant and would certainly warrant the derivation of more gen- 


eral differential equations. 
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DARIVATION OF PLATE COMPATIBILITY EQUATIONS 


1. General Discussion 


Borg, Hoppe and Kopchinski derived non linear strain 


compatibility equations in reference [4] by: 


a) First, assuming 63.= 0 
b) Second, neglecting entirely the requirements on 6z 
and simply utilized the non linear relations for 6x, by and 


Cxy as obtained from the non linear stress-strain relation. 
mj T=AI,E; t2Gn + ({1/-2mI.)E; +2mhr + sae 


The major limitation of the derivation in reference [ 4] is that 
it uses a stress-strain relation which neglects compressibility 
terms and other terms of the same order of magnitude as those 

retained. The more conplete stress-strain relation, as derived 


by Murnaghan (see appendix B) is: 
2 
(7) T=H=AL ili ‘F ((4-y I, 7 2mL.]E, +2(m ta-G)Lin tncoy + $607 


In addition, the derivation given in reference |4| assumes SY =e) 
as discussed on pages 13 and 24 which in effect says that the 
deformation distribution across the thickness of the plate must 
be constant and ignores any Poisson's effect across the thick- 


ness. This is not an unreasonable assumption for very thin plates, 


Sih 





but is a constraint that is open to some question for plates 
of finite thickness. (See the discussion of the order of mag- 


nitude of a on page 24 following equation (19) ). 


Using an approach similar to that used in reference cal 
Borg derived improved compatibility equations, the results of 
which are reported in reference [6]. In the improved derivation 
Borg used the second form of the stress-strain relation that 
includes compressibility terms of the correct order of magnitude. 
Borg's derivation reported in reference [6] is verified in the 
derivation shown later in pages 31 to 37 : It should be 


noted, however, that Borg's improved derivation assumes ow = O 
2 


and thus ignores any Poisson's effect across the tnickness of 
the plate. 

Intuitive reasoning would indicate that although the 
large deflection plate problem is neither truly plane stress 
(as implied by the von Karman theory), nor plane strain, with 
the exception of the very long narrow laterally loaded plate, 
it more closely approaches a condition of plane stress than 


plane strain. In linearized theory, a plane stress assumption 


would require that Y¥xz = Yy2 = 63 Om @,= Bury (C34 ey) #0 


Thus it would appear that a more valid derivation in the general 
case, for plates of finite thickness, would be based on au == (0) 


as shown in the plate strain tensor derived on pages 171 to 13. 
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2. Derivation based on generalized plate strain tensor 


ith Ow 
wl ae =o O 


By setting 62 = 0 and solving for qu in terms of 
W, X, Y, u, and v, we can obtain a relationship which is in 
accordance with classical large deflection theory and which 
will be valid in the majority of cases (one exception might 
be relatively thick shell pressure vessels subject to very 


high external hydrostatic loads, i.e. very deep diving sub- 


marine pressure hulls). 


Thus, solving for the terms of the stress-strain 


tensor equation 
a -G)L.n+nco +4+G 
iq te whe, + 26 +|(C-a)T, 2mI JE; +e 6) "( { (f 


2. 
we have (retaining terms up to the order € ) 
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sal = Cofactor matrix of the strain tensor 
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Retaining terms to the order of magnitude € , we have: 


~ /2w\* ow Ow O 
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|} gw ow -()° O 
ox 
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- qu Ya hs ft Fa (3 va) 
a ( ~ Da Naa 23 (ED (i hs ei 
yo es DPM tes ef 1 VSS ae 


2. 
i and Ly contain terms of order higher than 


and are neglected. 
Thus, writing out the stress-strain tensor equation, we have 
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The equations for the elements of the stress tensor can now be 


written 


vo on (eee HGNC] + 2 6 [8844 0) 
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Assuming 6x So as stated before, we have 


ba = 0 = (nv26) Me + nfde ade feat au + (z+ ote] 


solving for ow. 


Bey 

09) we (A du, av) em £252) (aus uy 
SZ At2Ge /\ Ox Sy 2(A +2G) 8 oy 

Mea peeceere CY Unich POR Steel 


HIG hy 


where Y = .3 is about .43. 


Thus ~ contains terms of the same order of magnitude as 
z 
the terms retained in the basic derivation, (€ ), and is of 


the same order of magnitude as the remaining terms in the first 
invariant of the strain tensor, and should, therefore, not be 
neglected. 

Substituting the expression for gy into the expression 


for by and by » we have 


gon) 6, = af Gale )-2285 (CONG) 
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Rearranging terms, we have 
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we can write the above expressions for 6, and 64 in terms 


of the engineering constants 
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Since we want to ultimately arrive at an equation in the form 
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as developed in the classical theory, we perform the indicated 


Gifierentiatvron. 
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deflection tensor. 
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the third term goes to zero and we are left with 
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Carrying out the differentiation in the last two terms, we have 
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If we consider the second invariant of the curvature tensor 
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and considering the expanded form of the 1st invariant of the 
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we can write the last two terms aS — 2c a Fe) |B, Gye + 2D Bhp 
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Thus, the equation for V F can be written as: 





eave = AY 3s +3) + (B+D) Sis as 2( D-c) one 


In order to get rid of the term A a +w) we go back to 


=o 


the stress relations 
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Now taking the Laplacian of the above expression, we have 
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Substituting this expression into 
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we can write the above equation in terms of the elastic constants. 
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3. Derivation of Plate Compatibility Equations based on the 


: QW _ 
assumption ae GO 
. re) ° . 
For the assumption 5>= 0 , we have the identical 


expression for the strain tensor as derived by Borg, Hoppe 


and Kopchinski [4]; 
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Thus, the terms of the stress-strain tensor equation 
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Terms ly and ia contain terms of the order a and can be 
neglected. 


Substituting the above expressions into the stress- 


strain relation gives 
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The equations for the elements of the stress tensor (equation 


[ 42 | ) are as follows: 
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Since, regardless of assumptions, we want to eventually have an 
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equation in terms of a stress function F or V F_ where: 
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we can carry out the indicated differentiation and add the 


resulting expressions 


Thus, we have 
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and the definition of the lst invariant of the large deflection 
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From the second invariant of the thin plate curvature tensor 
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by adding and subtracting -2 2 ow)" we have 
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Substituting this expression into equution (53) above, we have 
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(This form agrees with that derived by Borg [6] 


If we now make no assumption on 62 and work entirely with the 
meements of the stress tensor 6, , by and oxy » as was 


done in reference |4], we have, by adding equations (43) and (44) 
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Substituting equation (58) into equation (53), we have 
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(320 , No restriction on 62 ) 
Equation (62) is identical to that derived by Borg [6]. 
Thus, we now have three new compatibility equations under the 
assumptions indicated that can be tailored to a particular 


problem. 
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And finally, for comparison, the classical von Karman equation 
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DERIVATION OF EQUILIBRIUM EQUATIONS 


In the classical linear theory of elasticity, the equations 


of equilibrium are 








div [ +pF = 0 

or 
(64) ae + UT 4 Fe = 0 
(65) ys + 2024 + Py F, = 0 
(66) 9 Tis + SZe + O68 ae eee 


However, for the large deformation case, one must remember that 


the above equations are referred to the original coordinate system 


and that the curvilinear coordinates a, b, c, are identified with 


X, y, and z~ so that when the forces are projected, the changes in 
position of the points of the body due to its deformation are 
neglected. Stated another way, in the linear theory, no distinc- 
tion Vs made between undeformed and deformed values of the magni- 
tudes and positions of the elemental areas on which the stresses 
in projecting the forces, the rotation which 


act. In other words, 


an element of volume experiences as 
neglected. This assumption is open 
admissible in the general case. In 
ssary to take into account the fact 
c (the coordinates 


respect to a, b, 


tion). 


eee 


a result of deformation is 
to question and is far from 
the general case, it 1S nece- 
that differentiation is with 


of the points before deforma- 


Since the nature of the simplified linear equations of 





eguilibrium rests on the assumption of small elongations, shears 
and angles of rotation, and in view of the interaction of these 
quantities in the higher order equations, whether or not the 
non-linear terms can be neglected depends not only on the mag- 
nitude of the terms of the strain tensor, but also on the com— 
paritive magnitude of the corresponding terms of the stress ten- 
sor. It follows that the smallness of the angles of rotation in 
comparison to unity is not a sufficient condition for lineariza- 
tion of the equations of equilibrium. It is also essential to 
know whether the stresses which are multiplied by rotations are 
large in comparison with those stresses which enter linearly into 
the equations. The problems of elastic stability or of thin plates 
with large deformation are cases in point. 


The curvilinear form of the equations of equilibrium 
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can be transformed to the cartesian coordinates of the points 
of the body before its deformation. Novozhilov [2] has derived 


a vector form of the generalized equations of equilibrium as: 
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Resolving the vectors 6n, On, , bn, in the directions 1; jJ and lg 
we can write Novozhilov's equation 
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Projecting this vector relation on the x, y, z axes we have: 
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The above values of normal and shear stresses are not, strictly 
speaking, stresses. They can be called stresses referred to the 
dimensions of an element of volume before, not after, the deforma- 
fon. | 2 | 

Writing the non linear form of the general equation of 
equilibrium in terms of u, v and w (the displacements) we 


have: 
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Tue above relations clearly show the interdependence between 

the rotation and the stresses in the equations of equilibrium 

in the general case. However, the above equations are in gen- 

eral far too complicated to be of significant practical value. 
Since the basic assumption of this thesis is that the 

large deflection plate problem is both physically and geomet- 

rically non linear, the classical equations of equilibrium can't 

be used without a further analysis since their derivation is 


based upon 


(1) The assumption of plane stress 
(2) Geometric non linearity without regard to the physical 


non linearity. 


The classical equilibrium equations (appendix C) are valid 
approximations of the more general equations (75), (76) and (77) 
and thus, the geometrical portion of their derivation remains 


valid for the purposes of this thesis. Thus, the equation 


O* Mx ue 9° Mx 9’ M ere ow nA ow 
(78) Sn 2 2 ony 8 + 34 2} = (4 + Nx Sm WAN Ky Say Np > y= a 


remains valid, and the moments per unit length and twisting moments 


still are given by 5 
a 


(79) My, = a Z dz 


~$. ‘, 
aie My iia =, ann 


(831) May = Sf ‘Taq Bde 
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At this point, however, y depeavane from the classical 
approach becomes necessary. Instead of using the linearized 
stress relations to obtain a relationship between moments and 
deflection w, we will use the non linear form developed in 


mwendix B, 1.e. 


_ - 2 2 
(7) T=al,E, + 26n+|(t-ayL- 2m] Je, +2(m+a-e)ly + nce +46 
Using the analogy between the stress-strain tensor equation 

and the moment-curvature tensor equation as suggested by Borg 15 | 


we have: 


(82) M=AdcE, +26R +[ (L-a)de “2m. JE, +2(M4A-IM + NO +4+6R? 


where M is the moment per unit length tensor 


As “Mab -Q,c 
M = Mba Mob -QpC 
Ce 

- Qac -Qbc ‘a 2 


and (R is the curvature tensor given by ae > Aw 
ot i et 2 (SY) 113.68) on ) Re GS) 


(83) Y= fe aie ft = ARS (RY) +9. (ay) 2 (& io @w)+3.(2¥) | 


Yon = Mob Ve i 
fe te “e)  YRGDRGE) ARG BE]  &(52) 
4 ow a ow wow | owow 
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If we now use the classical large deflection theory assumption 
that the slope is small and the square is negligible compared 
to unity, and assuming Eulerian Coordinates, we have: 


_ | as ow tw ow 





E Ta re Sxt Ody xd 
(84) W%=2{ 1 -t 4 [= | gw Sw tw 
4x Ty Yue Sydx > uy 3) yS2 
sas gw Sw SW 
=r 2B a / OZ29x O2d4 92" 


At this point, it is necessary to consider what actually happens 
to the plate when we make the assumption Qu A O . This assump— 
tion, in effect, assumes a deformation gradient through the thick- 
ness of the plate (Poisson's effect) which alters the limits of 


integration in the expression for the moments per unit length. 


Thus, to be exact, for the assumption gw xO we should write: 
t, +5 
(85)  M, = [ oxzda 
-§-5 
45 
(86) { 
-t-§ 
©, +83 
(87) Mxy = Txy Zdz 
“ts 


where § is the change of the neutral surface due to the deform- 
ation gradient a across the thickness of the plate as shown 
2 


in figure 2 below. 
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FIGURE 2 


From figure 2, it can be seen that the total thickness is 
essentially not altered by either the assumption 3 > O or QW 0.4 
but that the plate deflection W _ is altered in the case of gv =0 
by an amount 4 , which in effect, requires an additional constrain- 
ing force F(z) such that the total deflection is now w= W(x, 4) +w(6)- 
Thus § is in effect, a measure of the error of the assumption 
ww = 0 
z 

A complete investigation of the effect of the § factor 
is beyond the scope of this thesis, however, it can be seen from 
equation (84) for the curvature tensor, that the curvature elements 
containing an are all of the order of magnitude et as dis- 
cussed on page 13. Accordingly, to the degree of accuracy of 


the other equations of this thesis, we can write the curvature 


tensor as t 





= 
2 a*w Y es —— 
ow Dx dy Y % Vey xt 
(88) = ww gi = ale Ay de. 
A Dyox y* o Vyx ry Yyz 
O O O a ali ele 
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From equation (88), we can see that we are left with the express-— 
ions for curvature that are identical to those contained in the 
classical theory in two dimensional form, where the third dimen- 
Sion terms are assumed to be zero. 


From the curvature tensor [equation (88)]| we have the 


following relations 





fa 


4 2 
fo) gic = -(> +o 4 +) Sly ap es 
ae a 
me 8. - tt fly wae invelvinguees 
de x fy (Z,) + terms involving Ye aie 
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2. 
= Ow Vy, aa 
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~) dy* Ox dy 


2 
(91) Sic contains terms of higher order than desired and thus 


can be neglected 





O O O 
(92) Cx k = 
O O Oo 
sw dw —(gw\" 
: © xt Oe Sn34) 


4 
= Terms of higher order than desired (€ ) 


ei 





4 
RK = terms of the order € and can be neglected 


4 
Ac K terms of the order €& and can be neglected 


Thus, we have, retaining terms to the order ¢€ 





(93) = saad ow A Mxo 
M = A ant Dys oO 1 0 + 2G 32w tw” 0 
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It should be noted at this point that retaining only those terms 
to the order Ee reduces the Moment—Curvature tensor to a form 
that is analogous to the linearized Hooke's law |= Waies +27 : 
Thus, we can see that the classical assumptions in effect force 
the equilibrium equations into a case of geometric non linearity 
and physical linearity where the third order elastic constants do 
not appear. This fact is born out by the fact that the third 
order constants 1, m and n are multipliers of non linear terms 
and cannot appear under the classical assumptions for the equi- 
librium equation. 

Inasmuch as we are now left with a linearized stress-strain 
relation and its analogous Moment-Curvature relation, the classical 
equilibrium equation derived on the basis of plane-stress becomes 
valid for cases where the actual case approximates plane-stress. 
Thus we have, from appendix C and reference 118] (neglecting body 


forces) 


4. -« PB, t/SEstw , af tw _ opSF Dw 
(94) Vow D a D &: ox? 7 Oxi 2 4* Mxdy Andy 


(PLANE STRESS) 
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Proceeding with the development for the assumption SW =o 
we have, from e= = in the two dimensional case and the 


analogy between tensor equations 


S) 
— 

(95) * i: 29/54 “ rSWe, —262 2 
Cox by oye —— 


Thus 
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(96) 6x = -z | a(&™ Sat) + 26 2% 
Be by = ~2 | a(S + sp) b2G So ye 
(98) Cry = ~Z 26 aay 


Substituting equations (96), (97) and (98) into the expressions 


for moments per ey: a we have 


sw 4 ot 
(99) mM, . [nz z*] a (xu+ =) +26 &™ | dz 


(100) My = JS ‘fa(z ees rr +26 2% Sy2] dz 





- M xy a is a 2g Sw dz 
ae KOY 
Performing the above integration, we now have: 
3 2 S aw 
no) My = —£ [A(Sa 138) 1268 J 
ere Se Sow 
(103) My = -=] a(3X+ Sy: +26 <% sat 
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Differentiating equations (102), (103) and (104) to obtain 


otMy _ 2 o Mxy 97M 
ra = Dn94 - ‘= 


, we have 


W105 ) UM ae ou 4 = + 2G ee | 








7M tw ot w 
(106) sah i -% BENG 344 ) + 26 eee | 
5 OM —- otw | 
mt O7 } 2 Say! ao 12 G dx7dy4* 
adding, we have 
3) QtMx _ My ie ee = oa otw tw 4, ow 
(108) VM 2 T+ & (2+ 26)( S55 + 23.587 Sap 


= -£ (A +26) VW 


Using th lationshi} = eGy So ECE) 
Sse e relationship A+t2G a on Ge ery 
= 2 
and G = a4) ,» we have A+26G = aK 1|-Y) 
| DP 2 
3 
Since the flexural rigidity 1s D= Et 
12 (1-2) 
We have finally 
2 > ye 4 
(109) oO My — 2 Mxy ,oM Mu ee Dh V 
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The equation of equilibrium for the assumption QW =o 


Se 


now becomes (neglecting body forces) 


D(I-¥) BW N,owW Nyy OW ow 
i” Gin ay) Viw=Rt oy de “83.04 7 t3y2 


or from the stress function F defined as 


(11) Nx = 6, 2 SF 





t OYyr 
(612 ) N = -~ > 
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We have 


4_ (i-2y) p £ (\-2y) (SF OW _ 50°F ow | ¥F Ow 
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As a check on the above equation, since the basis was a linear- 

ized tensor equation, we should be able to obtain the identical 

expression using the linearized plane strain equations resulting 
from the linearized Hooke's law 

ime ALE, + as where for plane strain 


we have the elements @z= Yuz = Y¥,xz = 2 =0 


—_ aie 


z= 


Thus, we have the expressions 


(115) 6x = Alex + ey ) + 2G ey 
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Pemee, 1m the linear relations 


oe "Ww ee =-2 2 
Vey = 5 7 = C% 3xdy where i =e ev * 34 
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we have the expressions 


(118) Sx —-AZ oe +o re as 262 UW. 


o) x2 





-2(2%4 Wee 343) “262 575 


19) by 


(120) oxy = —-262 (2%) 


Equations (118), (119) and (120) are Adlieifeal ke equations (96), 
(97) and (98), thus we could have derived the equilibrium equation 
(114) directly by the plane strain linear relations. The point to 
be stressed in the derivation of equation (114) is that the classae 
cal large deflection assumptions, when coupled with the assumption 


oa = O leads to a plane strain equilibrium equation rather than 


: . / 
a plane stress equation as derived by von Karman. 


In summary, within the limitations of the assumptions of 


classical large deflection plate problems, we now have two equilibrium 
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equations, the von Karman equation where the derivation is based 
on the linearized Hooke's law with strain a function of stress 

and the equation derived in this thesis with stress a function of 
strain. Thus we have an upper and lower bound on the assumption 


of plane stress or plane strain as follows: 


PLANE STRESS (Neglecting Body Forces) 








SS => 
——= = 


D Sy ox? a ax? o4™ OXd4 Bxdy 





(4) viws £ +e (REM, Brew pee Bw) 


PLANE STRAIN (Neglecting Body Forces) 


(114) Tiwe Uo2” Pt (1-2) (SF YW, SF OW _ FSF Otw 
(i-y)> Dp ~ D (i-y) Ge Sx2 — 9X* dy? Oxdy Sxdy 

As previously discussed, the actual case will be somewhere between 
the limits of the above equations. 

Applying the correction faye as a linear correction to 
the flexural rigidity in the data reported in references [14] 
and [19] and comparing the results with experimental results 
shown in reference [19] we can see by figure 3 that the correction 
appears to be about the correct order of magnitude and in the 
correct direction. Thus, for this case, it would appear that the 
actual case is closer to plane strain than to plane stress. On 
the other hand, one must be cautious about drawing hasty conclu- 
sions based upon figure 3 since there was insufficient data available 


in reference [19] to accurately evaluate the boundary conditions 


of the experimental results. Further discussion of the implications 


ye 








of the correction of Levy's solution and the divergence of both 


theoretical solutions from the Galcit experimental results is 


contained on page 62. 
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CONCLUSIONS 


In addition to the classical large deflection thin plate 


equations derived by von Karman, we now have three new compat-— 


ibility equations and one new equilibrium equation based on the 


assumptions indicated below: 


fs 


COMPATIBILITY EQUATIONS 

Classical von Karman Compatibility Equation [11] 

(2) VF = -Ed, = (y+!) 26 4: 

Based on 6% = 0, linear stress~strain relation and plane 


stress 


Borg Conpatibility Equations | 6 | 
4 
(>) V*F = -(G+8) duo - (4642) Sc 


Based on 62 = OF se = O anda non linear stress-strain 
relation 

ipl fm C= 2 al ale n 
(—) vee = | 21-Y) 4 IY Sirs j-v (4G + 2) Bac 
Based on SY = O, no restriction on 6z and a non linear 


SLress—Strain relatvron. 


Compatibility Equation Developed in this Thesis 

4 -~ mre = ie 
(d) YVRF= | 2) z = 26r] J... (7+1)(AG + 2) Gace 
Based on 6z = O, no restriction on ®9w and a non linear 


oZ 
Stress—-Strainervelatiton, 


=57= 





B. EQUILIBRIUM EQUATIONS 


Classical von Karman Kquilibrium Equation [11] 
4 E( SEN SFotw _ 4otF Siw 
(e) Vw = D ~ Sy ax? x2 Oye Oxdy DXOY 


Based on a linearized plane stress assumption (Assumes 





strain is a function of stress). 


Equilibrium Kquation Developed in this Thesis 
4, _ Ir2y p | t (-ay) (SF w , BE OW _ 9 VE Sw 
(f) VW = ap=p + dD deyyt ( dg2 ox * Sx dg? 3x84 3x04 


Based on non-linear development approaching plane strain, 
but linearized by use of the classical large deflection 


assumptions. (Assumes atress is a function of strain). 


From the above four compatibility equations and two equilibrium 
equations, we can now tailor the equation to suit a particular 
large deflection plate problem dependent upon the degree of 
accuracy required and the degree and type of deflection or 
loading involved. In other words, the equations can be tailored 
as follows: 

(1) Equation (a), or the classical compatibility 
equation can be used where the problem is linear physically, 
but non linear geometrically (moderately large deflections), 
or where the degree of accuracy is not critical. 

(2) Equation (b) can be used for very thin plates where 
the problem is non linear both physically and geometrically, 
but where the Poisson's effect across the thickness can be 
ignored, and the assumption 6z= 0 is valid. 


(3) Equation (c) can be used for thin plates where the 
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problem is non-linear both physically and geometrically, but 
where the Poisson's effect across the thickness can be ignored, 
and where the normal ($2 ) stresses can no longer be neglected. 

(4) Equation (d) can be used for moderately thick plate 
problems non linear both physically and geometrically, and 
where the Poisson's effect across the thickness of the plate 
must be included, but where the assumption 6z2 = 0 still remains 
valid (Thicker than Thin Plate Theory). 

(5) Equation (e), or the classical large deflection 
equilibrium equation, can be used in cases where the boundary 
conditions and loading are such that the actual case is closer 
to plane stress than plane strain. 

(6) Equilibrium equation (f) can be used for problems 
where a plane stress assumption is not valid and where the 
Poisson's effect across the thickness of the plate can be 
ignored (Sy = 0). 

Reflecting on the above equations, the question arises: 
why do the von Karman equations predict results that often closely 
approximate experimental results, but occasionally are signifi- 
cantly at variance? Particularly interesting is the fact that 
for the preliminary solutions Borg has achieved with equation 
(b) for circular plates with circular symmetry, the results are 
very close to those obtained with the classical von Karman equations 
(a) and (e) [7]. In view of the apparent wide variation in equa 


tions (b) and (a) 2 one musteasn y= eye 
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If we neglect the term in equation (d) involving the 
large deflection term Oup » we can note the similarity to the 


4 
classical von Karman equation i.e. 


(a) ae Mee) oe Sb. 
(dq) VF = — (+1) (4G +2) 326 


modified 


From equation (a) and (d) we can see that if the 


modified’ 
third order elastic constant n is equal to -4G, the von Karman 
equation is identical to equation (d) if the term [ mca) “2 -26| io 
can be neglected. 

As shown in Table I, Borg obtained values for n (corres- 


6 6 


ponding to Young's Modulus of 30 x 10 oe eo a) 1, 


Also, Smith obtained values for Austenitic Steel of ~40.0 x 10° 


and for 0.6 Carbon Steel of ~67 x 10° . kf we assume an average 
of the above values of n = ~47 x 10° and the commonly accepted 
meme Of G = 11.5 x 10° » OT 4G = 140m 10° » we can see that 


the actual numerical value of n_ probably does in fact lie 


very close to ~ 4G. Hence, we can approximate equation (d) as: 


VIF = [(2 +6) 1-2 Sao - (v4) 26 See 


The above equation still does not explain why the von Karman 
equation (a) provides valid results in the majority of the cases 
since we have only explained the right term and not the term 


involving Guo . 


If we analyze the function Sige = (2.4%, ies a) * +24)" | 


by vector algebra, we can see that grad Ww = oe S84 = Vw 
_ — 4 a 2 
. ~(QW) How)” - — 
and hence grad w- grad w (2¥) 43u) lv wl 
but the maximum directional derivative Ow. = | grad W| 
AS max 


= — 
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i 


ale 





cl ie | 


TAX 


hence, the vector meaning of Finp is that it refers to the 


Laplacian of the square of the maximum directional derivative, 
ee. vi | dw \" 
AS moy 

But in the classical large deflection assumptions, the square 
of the slope was considered small compared to unity, therefore 
the Laplacian of the square of the directional derivative will 
be much smaller than unity. Thus, within the limits of the 
classical assumptions for large deflection plate problems, the 
function involving Guo will be negligible unless the slope, and 
hence, the deflection is very large. Stated another way, when 
the deflections become such that the term involving Sino cannot 
be ignored, the classical assumption, in which the square of the 
slope is much smaller than one, 1s questionable. 

If we further consider the elastic constants associated 
with 4io ie Gee: -% - 26] and assume that n = -4G , 
we have Vane = ( ¥+6)(1-29) CA ~ (V+) 2G aie 


From Table I, for Borg's and Smith's values of m, we can assume 
an average value for steel of approximately m = -60 x 10° Or 


about —-5.2G. With this assumption, we can approximate the com- 


patibility equation (d) as — (166)0-2Y) Sup ~(Y+]) 16 Boe 


Soe 





we tory = 0,3 5 V*F = —.6446Gnp —2.66 952 or 


i = ~G [ 164 Sup + 2.6 Jac | 


Hence, we can see that the elastic constant contribution to 
the terms of the classical equation is about 4 times as great 
as the contribution to the term involving Sno . 

From Figure 3, it can be seen that experimental values 
for the simply supported square plate tend to diverge from 


theoretical values above ne of about 1.5 for the data presented. 


E 

Thus, from the hypothesis regarding the contribution of the term 
involving ic (eae and the results shown in figure 3, we can further 
hypothesize that for very large deflections (“s greater than 
about 1.5) the term involving Gio begins to become significant 
and the von Karman equations lead to increasing error with in- 
creasing deflection. Below the divergence area, the Ginn term 
is not significant and the von Karman equation is very close to 
the more exact equation (d). 

The classical equation of equilibrium (e) is identical 


in form to the equilibrium equation derived in this thesis (f), 


with the exception of the coefficient involving Poisson's ratio 


(1-2) | oe , 
Cre) ~ - For a Poisson's ratio of .3 the coefficient is 0.816 
and for Y = .25 , the coefficient is 0.89. However, as the 


deflections become very large and the deformations approach the 
plastic range, we have a different situation. In classical 


theory, when an isotropic material reaches plastic deformation, 
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Poisson's ratio is taken to be 4 |23]. However the coefficient 
OE goes to zero at Y = } and hence, equation (f) becomes 
Vw = O. Reference [25] states that some investigators have 
obtained values of Poisson's ratio in the plastic range much 
greater than } (as high as 0.8) which leads to further questions 
since the equilibrium equation (f) poeones We = -f (F,w,x,y). 
As a result, we are led to the tentative conclusion that equation 
(f) is only valid in the elastic range where the range of Y is 
0.25 to 0.35, or that for deflections approaching the plastic 
range, we can no longer ignore the Poisson's effect across the 
thickness of the plate (the assumption ox = 0 is not valid) 
and a three dimensional development similar to the development 
of equation (f) with sx # O will be required. It should be 
noted that the three dimensional development of an equilibrium 
equation where SY #z O , as proposed, would be based on the 
assumption that stress is a function of strain as opposed to the 
classical assumption of strain is a function of stress. 


One final point to be noted in the discussion of this 


thesis is that in equation (d) 


(a) vir = [2 G-2¥) “2 -26y | Sip - Cyt 4er2 )g,. 


the derivation was based upon the classical large deflection 
plate problem assumptions, but including terms applicable to 
plates of finite thickness, and with finite deflection theory. 
As a result, the development involves a compromise between thick 
and thin plate theory. Hence, we have a third category of plate 


theory which may be called "Thicker than Thin Plate Theory". 
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RECOMMENDATIONS 


The work done so far in the application of finite 
deformation theory to structures has hardly scratched the 
surface of a wide and virgin field of science. To attempt 
a discussion of areas requiring further investigation would 
be the subject of a thesis in itself. 

A few of the more obvious areas of interest to the 
structural engineer which require further investigation 
are: 

(1) To obtain more reliable values of the third order 
elastic constants 1, m, and n. 

(2) To obtain solutions to large deflection plate 
problems by use of the equations developed in this thesis, 
and for wiich there are experimental results available by 
Which one can obtain an idea of the significance of the third 
order constants, and to determine the validity of the hypothesis 
made in the conclusions of this thesis. 

(3) Determine values of § as applied to the limits 
of integration of the moment per unit length equations in order 
to determine the magnitude of the error involved in neglecting 
Poisson's effect across the thickness of the plate. 

(4) Through dimensional analysis, attempt to determine 
a physical relationship for the third order elastic constants 
as has been done with Lame's constants in linear theory. In 
other words, attempt to obtain a relationship between 1, m, and 


n, and the defining parameters of a material that determine its 
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thermodynamic solid state point (Temperature, Pressure, Specific 
Volume or Density, Conductivity, etc.). The long term impli- 
cation of the above analysis is the development of a combination 
Solid State Physics/Thermodynamics approach in the determination 

of material properties. <A logical extension would be to include 
time dependent functions for purposes of analysis of creep, fatigue, 
etc. 

(5) Attempt to determine equations for critical stresses 
for buckling or instability using the plate strain tensor derived 
in this thesis with a linear stress-strain relation in addition 
to buckling equations using the compatibility equations of this 
thesis which include the effect of the third order elastic con- 
stants of the non linear stress-strain relation. 

(6) Derive a Compatibility Equation under which no re- 
striction is made on either 62 or ow ( 6z £0 and 2 ZO ) 

oe Oz 
for use in "Thicker than Thin Plate Theory" where 6, must not 


be neglected (as in the casemou very deep-diving submarines), 
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APPENDIX A 


Strain Tensor i 


In tensor notation, the state of strain of a body, initially 
unstressed where point "a", having coordinates (a,b,c) , then 
strained to point "x" having coordinates ee) is given by: [5] 


dx* dx = da™ da. 


1 ) STATE OF STRAIN = 


Since dx - Jda, where J = O(x, 4, 2) ae Gbabanie ma brix of ee) 
a Ca, b,<) 
Ox 


Ox 

or Je 0a ob 

oY 

s Ss 3 
OZ 
ob 


oc = a tensor 


(ms 


00. 
* 
and d x = da’ J ie 


Ss is the transpose of the Jacobian 


where 
Sz 
2. Sense 
y*=/a l(x yz) =fax dy 
Ik 4 Sb db db 
5 x AF. 
Sz ae a4 Se 


and is also a tensor, we Can Say: 


STATE OF STRAIN = da” (I*T - E,) da 


(A-2) 

where 4 1s the 3 x 3 unit materi. 

Since X= GtK ; 5) = Denny , £ = C+w 

where * = X(G, bc) ; 4a 4 (a,b,c) and Z=2(4, b, < ) 





Ox _ Su 94 Iv z — ow | 
we have ee [ + ae , + = aoe : ge = | + oc | 
%¥ ~ QuQx - dU oy _ oV9 x say. 
$= 8% -3 ie Sy ce = (se , 
oz = Qu, dz = Qu 
oc ez 9oC¢ oc 
etc. (using chain rule for other terms) 


we Can express Jacobian as: 


a 
. = av OV IV 
(A-3) Re = aa | + Sb oc 
Ww gw. Peo 
5a ab aoc 
gu Vv oe) W 
+ 30 <i Se 
94 Vv Qw 
and x [+ = 
(A-4) Sp = 2c oe Se 
Q4 ov pes 
oc ac Yad 


hx panding a coal we have: | 
1 +224 4/34) gv)" ow Ju oe + 2Y OV, dude 4 Qwaw) [dy ,Sw OV Wy 24 du, Ow 
g " Si yz); +e] ~ Sate OG s 0Sp ob Soon) da. 9 aS 


(a-5) y*y= Fae ov ,avav, quay maw] 2 4 dv3V, dudy a 
3a dash) Oa a pOane ra $1425 eta yee) | [awades ise a 


ou dw 2 dvgv ,dudu ,2w dw 494 OY OV QU dU 4 ow dw fi+ oe te af] 
. sate song Nsaise porte ee Spe eaae” db d< * 3b | aoe 


a symmetric tensor 








Writing above expression in terms of orders of tensors, we have: 


is: — Es + a” + 2%, = E,+2(y,+y,2) + E3 $277 
where vf 1s a 1st order tensor 
and 2 is a 2nd order tensor 
(( is the combination 1st and 2nd order strain matrix 


and rearranging, we have 
(a0 fab 

(A-7) . You Yes 
| [<o Yel 


= = “en se 2 
ee 





From the expansion of J*J and the relation “) = 5 Cie E3) 


we have 
* (3+) 4 Qty) 
(A-8) "= 2 (3¢ +3) ay L(o¥ + oY 
(52 *3%) x (34 +35) Se 
(3eyegz) ay] Seepage gray] [Bebe er +gre] 
+ f[Seser Ree N] Be) Ge eGey] ese earae 332] 





The above expression expresses the strain tensor in the Lagrangian 
coordinate system where the movement of each particle in the body 
is followed (i.e. the initial coordinates of a particle were (a,b,c) 
then the final coordinates of the same particle in the deformed 
mosition at (x, y, z) is xvanbye pu eemc uc byc) |) 21>) 
If the assumption is made that the deformation is sufficiently 
large that the 2nd order strain tensor must be included, but is 
sufficiently small that "a" is essentially the same as "x", the 
expression for strain can then be expressed in the Eulerian coordi- 
aeane system (i.e. conditions are expressed at each point in the 
deformed body, and a particular point rather than a particular 
particle 1s considered and quantities considered are functions of 
noe and Zz ?). 


Hence, we have the approximation 


a ub > 
e x ‘ Yx4 £YVxz on 23% + $y ) 5 ($+ wy 
P= [a%e ey Yaa 134) 7(34 + $Y) 
' 
> Yzx Ly. ez I/qu td 
a ere) a) OF 


(using American notation for elongation and shear strain). 
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APPENDIX B 


Derivation of Finite Deformation Stress—-Strain Relation 
I. The fundamental relation of elasticity theory connecting 


Stress and Strain 2s: 


(B-1) T =(B) J oF i | [1] 


where T aus the stress tensor 





6 x Cx Cxa 
es d 7 OF = sacobian = 9(X, 2) 
fax by a ) 9 (4, b, cy 
Cex Coy ba 
(2% is the compression ratio = | 
det J 
Dp = P(N) = energy of deformation / unit 
initial vol. = a W 
/ = Strain matrix = 


( B-2 ) / = : (T*T - E3) 
(see Appendix A) 


Il. If the deformable medium is assumed to be isentropic, it 
is elastically insensitive to every rotation of the initial 
Cartesian reference frame and hence Py) = D ( Rv R) 
where R = the rotation matrix. 
Further, a deformable medium is elastically isentropic if and 


Only 1h P (1) is a function of the three invariants I 


I, of ( . 


1? I, and 





In reference |1] it is shown that: 





j Oo o 
(B= 3) oi = Es = O | O 
ob 0 Oo ] 
(B-4) ot ae ae 
( 
(B-5) ols = 
B~5 7 Co / 


It follows, that for an isotropic medium 


9g of] 9291, , SP OF 
6 = 3057 TSK 3p * Sz oy (chin Rule) 


or Og 
7) 


MI 


see, + 28 Les-) + oe 2 


Ill. Assume the energy of deformation per unit initial volume 
is developed as a power series in terms of the three strain 
invariants I, 1, and I... Jtak ” is assumed infinitesimal, 
I, is of the @rdervol Cael of the order Ee” ytd 


es 
I, of the order €° ° We can then say: 


3 
fe) Py) = PiCipee ee = Ae A,+AL+4A,L oT 
+ G,l. + 6,1, bree 
+C,T, +Cilpee + DUT 


3 
Neglecting terms of higher order than € we have 


(B-8) P(1) = A, +A,I, + LAL, + J Ae + BT, +C0l., +0 
FA 





or in terms of GB = Bg, + GD, + DM, t+ Py teh 
where QD, consists of terms independent of v4 
g, terms linear in elements of Y 
, terms quadratic in elements of 7 
and hence a linear combination of ite aa. and 


D, a linear combination of jae ip EN and 1; 


3 
x § oa L Ci 
9) P29 + At +(AVe81) (eb + ents + DT) 
P, Do D3 
Differentiating with respect to " » we have: 


(B-10) 9% =A,ol, 4 A,I.oL 4 8, 9h + AST, OL 





ola Ss ‘ ol. = l Ee = ‘ ols = 
Since oy” Ean. 3y ( 3 y } ) 7 a 
a a i 2 
(2-11) s = AE, +A, LE, + 8, (Le;-¥) + > TEs +& TEs +1,(6~4)} 
+ D, Coy 
Let 
_ ve 

(B12) (Aat+8&) = A = GER Lame's Constants (A&M ) 
(B-13) B, = -2Z&Ar = -26 ( Y = Poisson's ratio) 
(B-14) Ast = 

h 
(B-15) Co = —~2wym l,m, and n = 3rd order constants 
(B~16) D, = n 





Substituting the elastic constants into 2@ we have; 
4 22 = = : 
(B17) 7 = AEs +XLEs +26y (LL, -amTz) Ey tamIiy +n coy 
which is the 3rd order approximation obtained by Murnaghan Ene 
The higher order relation between stress and strain 


x 
resulting from the generalized relationship T = Se (x7) J 


is then: 
(B-18) T= i TI Ase + ALE, + 2ey +(LLM2mL JE s+ 2m Ly +n coy | dea 
If the elastic medium is initially unstrained, ox Ss j/ eats. fl 
is a zero or null matrix, therefore J is a rotation matrix only. 
If To is the initial stress in the unstrained condition, the 
above tensor equation becomes 
(B-19) To = AiE3 
where Ty is due to a hydrostatic pressure or tension. 

For an isotropic medium, we can say A, = sd hoe 
initial hydrostatic pressure) or TL = c-P,E, - With an 


initial hydrostatic pressure of af equation (B-18) becomes: 


(B-20) Ts os T|- poe, FALE,*Gy + (AI>-2 mL, Es ey + ncoy | hg 


For an isotropic medium, the coordinates of the stress tensor 
are furnished in the final rectangular Cartesian reference 


frame in which the coordinates of ape are furnished by the 


elements of Mz J *J by the elements of the matrix /% mM 9é 
Pa. Oo” 


Puere J =R Mie and T*J=M= 27+ 


Thus, the relation between stress and strain becomes: 


B-4 





oan) Te & (Essay) Pres HOLE, #2ey) #(An ames 


+ ae a9 rina 


T" -2 |-PoE s+ AT,E, + =) + (AL. am1, JE, 1 oes 


+ oe + bie | 


oe 10T +P o =O 


(B22) T= 2 [AL Es+2en4 (AL -2m1,)Es pan h? ina : 46n7 | 


Murnaghan [1|, has shown that 


(B-23) (oe det (Es +2y) = (det J) = }421,+41,+8!3 


OB & 


-, 
ee (tear +4I,+91,) * 


To a second order approximation, the compression ratio 1s 


(B-24) SE = (1421, Pai) Su-7 +@1,°-21, ) 


neglecting terms higher than 2nd order, we get 


T 2 AL, Es +2Gy + (AT,*-2mli Es + 2Cmeaj | ss ey +4607 — ALE; “261 y o 
(B25) tear Ey + dey + LCA-0* -2mT, | FE, + 2(mtaA-G) liy 


+n "l re ae 





SUMMARY 


Third-Order Approximation of stress-strain relation (three 


dimensions). 


(p-25) [T= ATE ¢264[(L-2)1-2m JE, + 2(m-+2-6) Ly +n 7 be 


where ale = stress tensor 
fl = straenvensonr 
E 3 = identity tensor 
rie = First envibiant of strain tensor 
es = second invariant of strain tensor 
A G& = second order Elastic Constants 


(elastic constants of Lame) 





fim, n =) Uhind Crder  ilastre Conus tants 
A = _YE _ 26? 
(i+ ¥)(1-2Y) |-2y 
| » 
= Ire : = 
4 oisson's ratio XatG) 
EE = Youns's Modulus = G(3At2G6) - 26C1+Y) 
At+G 
CS = = = Shear modulus 
2(1+¥) 





APPENDIX C 
DERIVATION OF CLASSICAL EQUILIBRIUM EQUATIONS 


IT. Balance of Forces 


Considering the elemental plate shown in figures C-1 and 
C-2, we have for the x and y component forces 
x component 


mC=1 ) (Se 2 dx Jd - Nx dy r (Nyx fl ae Ay) Ax ~ Nyxcke + KArdy =0 





= iu dx dy fe og +X dx dy o 
= SNx 4 INyA 4X = 0 
9 xX dy 


where »'@ is the x body or tangential force component per unit 
area of the middle plane of the plate 


y_ component 


Similarly we have for the y component forces 


(c-2) Ny , ONY yy so 


gw 


If we assume SEs and = are small compared to dx and ay , 
we have for the projection of the Nyx forces on the 2 axis: 
2 
(c=3) aN sw yow oY = ow 
CN + 9Na dx (3% Caan r)ay N x xe 


— 9Nx ow N, 2 AE 2 
= ae ou dxdy + & 54 olxdly + ON aw Aik dy 


CI 





or, reducing to second order terms, we have 


(C-4) Nx O*w cl x Ay 
9 xX* 


and similarly, the 2 component of the Ny forces becomes 


(C-5) 4 ow 
Saye aan 


since Nx = Nux from equilibrium of moments about the 2 axis 


the 2 component of the Nxy and Nyx forces is 


Oo) 2 Ny Sw 


Thus, the membrane portion of the 2 component forces become 


(C-7) Nx ay + 2Nn S¥ S134 Ba ‘one O*w 


S44 


The effect of shears Q, and Q, », aS shown in figure C-3 also 


contribute to the equilibrium of 2 component forces, and can be 


shown to be: oQx DQ 
— + aa) ) dx dey 





The lateral load component is P dx du and the body force com- 


7 t = ow nd — Qw 
ponents are X $Y Axdy a NOR _ 
Thus, the combination of 2 components in the a, case is: 


Po) Dt ae pe + Nx Ww + aN +N Se tw -xw “Y 5422 
‘ i 





BALANCE OF MOMENTS 


Considering figure C~3, by taking moments about the x axis 
and considering the right hand rule for positive moments, we 


have 
(C-9) (May + oMry di) y — Mig oly + My Ax _ (My + Stak) dx 
+ ( Qu 1 SU oly) Ax dy =o 


Simplifying and retaining terms up to the second order only, we 
have 
(C-10) O9Mx, — OMy 4 Q) 
> X © 5 
Similarly, by taking moments about the y axis, we have 


(c-11) 9Mx , OMyx - Q, =0 
Ox Dy 


Differentiating the 1st equation with respect to y and the 


2nd with respect to x and adding, we have: 


C-12) O°Mx 4 O*Myx 4 otMy — O'My Ke a 
or DuSy _ sye roy ra ax si 











Or Since My x = = Mxy (Oxy = Tyx) 
we have 
2 Mx a*M —~— O9Qx 4, 9A 
13) 2S. — = , Soy 
, ox? 3x4 r DKSY ox 20 
But from the 2 component force equation, we have 


(C-14) OQx_ = Ow Oy ow — Sas 
— Fra = -(p + Nx Tso," NG 94 KX QW Ye ) 


C—3 





Thus we have 


QM _ 9 OMay 4 My. _ (pp tNxd 4 ANS 4 Nyt —K ow _ Yow 
o-15) sx 2 Snay Sue = (p onras “Gxdy Jay? XS TS 


The moments resulting from the stresses distributed on the x and 


y faces of the differential element shown in figure C-3 are 


C/2 
oy = I Be ae 
-t), 
OC) 
ee=17 ) PO J b,zd2 
~o/, 


(C-18) M ig 
a2) =a Ge zZdz 
~ty J 

From the plane stress equations of the linearized Hooke's law, 
we have 

(C-19) @ = ( 6x ~¥ by ) 
(C-20) @, = + (6 eG) 

J En aay x 


e211) Yay > 2(1+yY) ny 
E 


or, in terms of stresses, we have 


fC-22) bx = = (€, + ve, ) 





since 








_ _ 2z 9?w 
Ce. = 2 _ Z ox? 
= _ = 3?w 
" 
Q*w 
Yru = ~ 22 dyay 


by substituting the stress relations (C-22), (C-23) and (C-24) 
in terms of Z and the curvatures into equations (C-16), (C~17) 


and (C-18) and integrating, we have: 





e—2 5 ) Mx = —D ve ta oH 
So” bs 
(C~26) My =D ( oy: * ew 
(C-27)  Mxy = D(i-y) S7W 
dxd4 
where D is the flexural rigidity dD = Et, 
ee 


Differentiating equations (C-25), (C-26), and (C-27), we have 


ed DL dx Say) 








iM . pn / SY -y oe 
(C-29) ay a D( Sys * Y Soy) 


(C-30) 9° May = DCt-y 
Sx Sy ) 2%) = 


¢ 





Thus 


31) S'Me _2¥ May, YMy = DOW pry dw + PW 4 2(rv) dee =-DVW 
(C-31) ae 2 Say Se oe +2 BOG Oy ot) 





C-5 





Finally, substituting equation (C-31) into equation (C-15) we have 


(=32) Vow = £ : (MB + any Bi, + Ny Sw -Xqw “TSy 
ay? 


From the Stress Function F defined as 








oF = by = INx 
Ou ~ 
6 
Se ae ~ Ny 
oO yr Ou Lz 
S22 = N xy 
Sx OY oxy i ss 


We have 





(233 } Vw = 4a. oi o°w a QoF Ww + SF 3 'w x 
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Figure Cel 


= Undeformed Plate Element 
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Figure C-2 
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